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^ ■ Abstract 

a' 

O ' We use microscopic many-body theory to analyze the problem of itinerant ferromagnetism in a 

O 

repulsive atomic Fermi gas of Hard Spheres. Using simple arguments, we show that the available 

theoretical predictions for the onset of the ferromagnetic transition predict a transition point at a 

CN ' density {kpa ~ 1) that is too large to be compatible with the universal low-density expansion of the 

CN ■ energy. We present new variational calculations for the hard-sphere Fermi gas, in the framework of 

rn . 

^^ , Fermi hypperneted chain theory, that shift the transition to higher densities {kpa ~ 1.8). Backflow 

(N 

correlations, which are mainly active in the unpolarized system, are essential for this shift. 
^ ■ PACS numbers: 67.85.-d, 05.30.Fk, 75.25.-j 



I. INTRODUCTION 

Experimental and theoretical progress over the last decade has allowed for a deeper 
understanding of phenomena such as superfluidity and pairing along the BEC-BCS crossover 
in degenerate Fermi gases. [l\ In these studies, the diluteness condition of the system allows 
for a characterization of the interaction by its s-wave scattering length a, a quantity that is 
easily tuned by making the system approach a Feshbach resonance. In most cases, however, 
the underlying interaction between atoms is short ranged and attractive, while little attention 
has been paid to the case of short range repulsive potentials. 

The role played by short-range repulsive interactions has been revived more recently, with 
a series of experiments discussing the possible onset of a phase transition to a ferromagnetic 
state in two-component ultracold Fermi gases. In an initial work, Jo et al. [21 describe an 
experiment with ^Li atoms which is supposed to implement the Stoner model 3|, \^ of mag- 
netic interactions, and find clear signatures of the ferromagnetic transition by monitoring the 
total and kinetic energy of the gas and its volume. They conclude from these measurements 
that a transition to the ferromagnetic state occurs around x = kpa Ri 1.9 ± 0.2 with kp the 
Fermi momentum of the unpolarized system. In a second and more recent experiment |5|, 
the existence of a ferromagnetic transition in the same system is ruled out, arguing that 
very fast molecule formation leads to local heating and to a pairing instability, destroying 
the ultracold nature of the gas. They even conclude that short range repulsive interactions 
with fermionic species can possibly not be realized in nature. 

The existence of a ferromagnetic phase in dilute, ultracold Fermi gases interacting through 
short-range repulsive forces has also been discussed in recent theoretical works, triggered by 
the previous experiments 2|, |5| . Positive scattering lengths can be achieved by moving along 
the upper branch of a Feshbach resonance associated to an attractive short-range interaction 
once a bound state is formed, or by increasing the range of an overall repulsive potential. 
These two mechanisms lead to different physical situations, as pointed out in Refs. p, l7|, 
when a description based only on the scattering length is not enough. However, in both 
cases it can be argued that when the system is so dilute that the scattering length and 
the range of the potential are small compared to the interparticle distance, only s-wave 
scattering processes are relevant. Under these conditions, all higher order partial waves can 
be discarded in the description of the many-body syste. In any case and motivated by the 



previous experiments, the question of whether a ferromagnetic transition takes place in ul- 
tracold Fermi gases has achieved revived interest recently. However, no consensus has been 
reached yet. For instance, in Ref. [Sl-ll 1| the authors find that a transition exists in different 



systems of harmonically trapped repulsive fermions, while similar conclusions regarding in- 



finite systems are reported in Refs. la, 02, 13 1 . In other works, the existence of a transition 
is either ruled out or argued to depend on the details of the system analyzed [1^, 115| . In 
this work, we address this problem again and analyze the possible onset of a ferromagnetic 
transition in a system of fermions interacting through purely repulsive, short ranged forces. 
We believe that the question of weather a system of repulsive fermions can undergo a fer- 
romagnetic transition is still relevant even though in actual experiments the high x limit is 
usually reached by moving along the upper branch of a Feshbach resonance. 

II. MODELS 

In the case of a fully polarized Fermi system, s-wave scattering is forbidden by the Pauli 
principle, and thus the total energy equals the corresponding kinetic energy of the underlying 
Fermi sea if the density is low enough. In a non-polarized medium, s-wave scattering between 
particles with the same spin orientation is also suppressed but not between atoms of different 
spin. All that information can be collected in the following, effective opposite-spin-channel 
(OSC) Hamiltonian 

with m the mass of each of the A^ = A^-^ + A^^ atoms, A^ and A| being the number of 
particles with spin up and spin down, respectively. The Stoner model |3|, \^ corresponds to 
a mean-field approach to this Hamiltonian when the interaction is replaced by a two-body 
pseudopotential with a coupling constant proportional to the s-wave scattering length a. 
As commented above, this model has been taken as the starting point in recent theoretical 
analysis of the ferromagnetic transition. Contact interactions in standard Monte Carlo 
simulations of dilute quantum gases are usually replaced by model potentials of finite range 
tuned to reproduce the desired values of scattering parameters. One of the fully-repulsive 
potential models commonly used to study low-density properties of degenerate quantum 



gases is the hard-sphere interaction, 

. s 1 +00 r < R 
Vir) = { - , (2) 

otherwise 

with an g-wave scattering length a equal to the core diameter R, which been used recently 
in Refs. [d, M, 



16|. 



When a finite range potential is explored, higher order partial waves produce contributions 
that show up when the diluteness condition of the gas is not fulfilled. In this case the OSC 
model, which explicitly restricts this possibility, shows its limitations. When applied to 
the experiment with ^Li atoms of Ref. |2|, all models predict a ferromagnetic transition 
at a point where the system is clearly not dilute. With kp the Fermi momentum of the 
paramagnetic phase, the Stoner model predicts the transition at kpa = 7r/2, while the 
Monte Carlo simulations of Refs. p, |l3| lower that prediction to kpa ^ 0.9, with minor 
differences between the case of a fully repulsive interaction and a short-ranged attractive 
one. With a transition point at kpa ~ 1, it is no longer evident that the OSC model can still 
be used to describe the system. In this work we study the possible onset of a ferromagnetic 
transition in a spin 1/2 Fermi gas of Hard Spheres (i.e. a two-component Fermi gas) using 
existing low-density expansions of the energy per particle and a variational approach in the 
framework of the Fermi Hypernetted Chain Equations (FHNC) [l8|,[l9|. As we do not restrict 
the analysis to the weakly interacting regime, we allow for the contribution of scattering in 
all partial waves by adopting a more general Hamiltonian of the form 

where indexes i,j,.. and i',j'.. label spin-up and spin-down particles, respectively. We 
stick to the case where the two-body interaction does not distinguish between different spin 
configurations, but explicitly include interactions between spins with the same orientation, 
in contrast to the OSC model. Since we are analyzing the gas of hard-sphere fermions, we 
use ([2]) for V{r) in all three channels (up-up, down-down and up-down). 



III. RESULTS 

We are mainly interested in the comparison between the energy per particle of the para- 
magnetic, unpolarized system [N^ = N^ = N/2) and the ferromagnetic, fully polarized one 
(A^-^ = A^). The density determines the Fermi momentum, which is different in each case due 
to the spin degeneracy. For the polarized and unpolarized systems one has kp = (Gvr^p)^/^ 
and kp^ = (Stt^p)-"^/^, respectively, so kp = 2^^^kp^ . As in previous works, we set kp = kp^ 
as the momentum scale unit and plot our results in terms of the dimensionless quantity 
X = kpa. Due to the spin degeneracy, the paramagnetic phase is preferred in the noninter- 
acting a — )■ limit at zero temperature. The stability of the system when the polarization 
A = (A^^ — Ni)/N increases is determined by the inverse magnetic susceptibility x~^j 

The critical density against spin fluctuations corresponds to the point where x~^{^c) = 0. 
On the other hand, a sufficient condition indicating that the ferromagnetic phase is preferred 
is given by the criteria E^{x'^) — E^^{x'^) = 0, with E the energy of the system. This second 
estimation, which is the one we use in the present work, is fulfilled at a slightly larger 
density than the first one (xc < x'J. For instance, in the Stoner model Xc = 7r/2 2^ 1.57 
while x'^ = (97r/10) (2^/3 - 1) - 1.66. 

Several expressions for the ground state energy per particle of a system of interacting 
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fermions at low densities have been derived in the past 2[H22l|. All of them rely on per- 
turbation theory in terms of a renormalized, effective interaction obtained from a G-matrix 
which contains information about the multiple scattering of particles moving on the corre- 
lated medium. For a strong repulsive potential the massive summation of selected terms is a 
convenient strategy to get a rapidly convergent perturbation series, and becomes absolutely 
mandatory for the extreme case of the hard-sphere potential due to the divergent character 
of the matrix elements of the bare potential. A very useful expansion containing terms 



beyond s-wave contributions is the one derived by Bishop in Ref. 22], which is expressed 
in terms of the s- and p-wave scattering lengths a and a^, respectively, and tq, the s— wave 
effective range, as 
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FIG. 1. (Color online) Energy per particle of the polarized and unpolarized hard-sphere Fermi 
gases. The black solid, violet dotted and green dot-dashed lines show the low-density expansion 
of Eq. ([5]) for the unpolarized system, the same with the p— wave contribution (proportional to 
Cp) removed, and the universal prediction of Eq. ([8|). The blue dashed line depicts the prediction 
for the polarized system given in Eq. ([7|). The red dot-dot-dashed straight line is the constant 
value assigned to the polarized system in the OSC model, while the red squares and solid black 
circles stand for the DMC calculations of Refs. [6j and [13l | for the unpolarized system in the same 
approximation. 



E 
N 



2m 
1 



+ 0^ 



2,, , 4(ll-21n2),, ,2 



+ 



IOtt 
1 



{kFro)ikFay + (0.076 + 0.057(z/ - 3)){kFay 



Stt 



u + l){kFap)^ 
16 



+(z/ - l)(z/ - 2)— -(47r - 33/2)(A;^a)Mn(A;i.a) + ... 
277r'^ 



(5) 



with z/ = 1 (2) the spin degeneracy of the polarized (unpolarized) system. In this expansion, 
the p-wave scattering length is defined in terms of the p-wave phase shift Si{k) as a^ = 
(— fc^cot(5i(A;)/3)^^/^| r,- For the particular case of the hard-sphere interaction, where a^ 



a and tq = 2a/3, the previous expression reduces to the one originally derived by Efimov 20|. 
For an unpolarized system of Hard Spheres and with x = kFCL, the previous expression 
reduces to 
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while for the polarized system one recovers 

— = — — I X in 

A^ 2m [5 Svr J ■ ^ ' 

Notice that in these expressions and up to order x^, E^^ /N contains s- and p-wave contribu- 
tions besides the kinetic energy of the corresponding free Fermi sea, while s-wave effects are 
suppressed in E^ /N. However, there is still a p-wave contribution of order x^ that becomes 
important when the density increases. This can be seen in Fig. [H where we compare the 



E 



FN-DMC results of Refs. 13|] (black circles) and [6| (red squares) for the unpolarized gas to 
different instances of Eq. ([6]) obtained by keeping all terms (black solid line), removing the 
p-wave contributions (violet dotted line), or keeping only the universal terms (green dot- 
dashed line), which are of order x^ and depend on the interaction only through the s— wave 
scattering length 

(8) 
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To this same order, E^ /N reduces to the kinetic energy of the underlaying polarized Fermi 
sea (red dot-dot-dashed line). 

The previous expansions are asymptotic in nature and so it is difficult to accurately de- 
termine its convergence radius. However, in the particular case of a hard-sphere interaction, 
they are known to reproduce the equation of state up to relatively large values of x. Fur- 
thermore, they can not be directly compared to the OSC model as it is implicitly assumed 
in their derivation that the Hamiltonian contains the same interatomic potential in all chan- 
nels. It is however possible to obtain a prediction for the OSC model if the leading terms in 
the expansion are renormalized accordingly. The OSC model discards interactions between 
particles with the same spin orientation, which in a first approximation can be thought as 
forming pairs in a triplet state of total spin S = 1. The wave function in configuration 
space for these pairs is therefore totally antisymmetric and thus s-wave scattering processes 
between them are suppressed. In this way and to leading order, the s-wave contributions in 
Eq. iQ should be the same in the two model Hamiltonians of Eqs. ([1]) and ([3]). For 5* = 1 
states, p-wave scattering is the leading contribution, but only the symmetric combination of 
unparallel spins is allowed in the OSC model. Assuming all three spin configurations con- 
tribute the same to the triplet state in the Hamiltonian H in Eq. ([3]), we conclude that the 
p-wave contribution to Hqsc should be one third of the p-wave contribution to H . Taking 
all these facts into account, we can build an expansion of the energy per particle of the OSC 



model by simply weighting with an extra factor 1/3 the term proportional to a^ in Eq. (Q. 
Doing so leads to the green dashed line in Fig. [H which accurately reproduces the Monte 
Carlo calculations of Refs. |6| and [l3| which were obtained from the model Hamiltonian 
Hqsc- Accordingly, one may conclude that the low density expansions given above are valid 
in the whole range of x values considered, and that both s- and p-wave scattering processes 
contribute significantly to the total energy per particle when x increases. 

Assuming then that the above low density expansions provide a good description of 
the energy per particle of the polarized and unpolarized systems at x ~ 1, one can draw 
several conclusions from the different curves in Fig. [H For instance, one sees from the figure 
that the x^ terms add important contributions to E^^ /N. Adding only x^ terms to the 
s-wave rises the energy per particle slightly, but the most significant effect appears when 
the p-wave contributions are added, bringing the violet dotted line into the black solid line 
corresponding to the prediction for the full model of Eq. ([3]). Therefore, adding p— wave 
effects is more important than dressing the s— wave contributions by including the effective 
range terms. Furthermore, the good agreement between the Monte Carlo predictions and 
the properly modified low-density expansion for the OSC model, indicates that while p-wave 
contributions are important, no higher order partial waves contribute significantly, as these 
are not included in the expansion. We thus conclude that s- and p-wave processes are the 
only ones that contribute significantly. 

It is also apparent from Fig. [1] that the universal behavior, where all predictions yield 
essentially the same, ceases to be valid already at xq ~ 0.4, quite below the value x ~ 1 
where the simulations based on the OSC model predict the transition to the ferromagnetic 
phase. The figure also indicates that the energy per particle of the polarized gas starts 
to deviate from the constant prediction, valid when only s-wave scattering is considered, 
at about the same value x = Xq. Once again, at x ~ 1 p-wave contributions introduce 
significant corrections that can not be neglected. We thus conclude that, at least for the 
hard-sphere Fermi gas, the low density expansions set the limit of validity of the OSC model 

at X ~ Xq. 

Knowing the limits of the OSC model, the next step in our discussion is to analyze the 
results provided by the full Hamiltonian H of Eq. ([2]) with the interaction acting in all the 



spin channels. To this end we perform a variational 



Fermi Hypernetted Chain equations (FHNC) [18|, 
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calculation in the framework of the 



23|. Although not exact, this method 



has been successfully used in the past to describe strongly interacting liquids like pure 
^He 24J-|26|. ^He-'^He mixtures |27| and nuclear matter 28|, |29|, and is therefore expected to 
accurately decribe the physics of the present problem, We use a variational Slater- Jastrow 
wave function of the form 
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where /(r) is a two-body correlation factor, while D^ and D^ are Slater determinants of spin- 
up and spin-down plane waves filling momentum states up to the Fermi level, and including 



backflow correlations 



30|. In our FHNC calculations, the sum of elementary diagrams is 
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25| . The results obtained 



approximated using the interpolating equation approximation 
with this wave function iQ are compared with those obtained from a less sophisticated 
version \1/ where backflow correlations are removed. In both cases, we take the Jastrow 
factor /(r) to be the same for all spin channels, and equal to the solution of the optimal 
Euler-Lagrange Hypernetted Chain equations (HNC/EL) (23l.l3l| for the underlying Bose gas 
of atoms of the same mass m interacting through the same potential and at the same density. 
We have checked that, by using the optimal f{r), the energy per particle improves noticeably 
compared with other, simpler forms containing few variational parameters, specially in the 



case of the unpolarized gas. On the other hand, backflow correlations enter exponential 
the Slater determinants of plane waves through the renormalized position coordinates 
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with A, To and ujq variational parameters that are optimized at each density. Backflow 
correlations are capital in an accurate description of the energetics of strongly correlated 
systems like ^He, 32|, |33| and have been reported to have also a non- negligible impact in the 
problem considered here at the densities where the ferromagnetic transition is predicted |6|. 
Figure |2] shows our results for the total energy per particle of the polarized and unpolar- 
ized systems in different approximations. The solid and dotted lines show the low density 
expansions of Eqs. ([6]) and ([7]), respectively. The open triangles and diamonds correspond 
to the predictions for the polarized and unpolarized systems using the HNC/EL-optimized 
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FIG. 2. (Color online) Energy per particle of the paramagnetic and ferromagnetic phases in different 
approximations. Solid and open diamonds: FHNC results for the mipolarized gas including or not 
backflow correlations, respectively. Solid and open triangles: same for the polarized system. The 
solid and dotted lines are the predictions of the expansions in Eqs. ^ and ([7]), respectively. 



wave function of Eq. ([9]) without backflow correlations, while the solid symbols stand for 
the same predictions including backflow. Several conclusions can be drawn from this flgure. 
On the one hand, it is remarkable the fact that while the low density expansion for the 
unpolarized system is everywhere close to the full many-body calculations, this is not the 
case for the polarized system, which overestimates the total energy per particle of the gas 
when X > 0.6. On the other, while the low density expansion predictions do not cross in 
the range of densities analyzed, our many-body calculations show that the energy per par- 
ticle of the polarized system decreases at a higher rate than the corresponding one for the 
unpolarized gas. It is also remarkable to notice that the total energy per particle of both 
systems noticeably depends on the quality of the wave function employed. As expected, the 
one including backflow correlations provides lower variational estimations. In any case, the 
inclusion of backflow correlations not only lowers the energies, but also pushes the ferromag- 
netic transition to a higher density, identifled by the point in the flgure where the energy of 
the paramagnetic state equals that of the ferromagnetic phase (a sufficient condition that 
determines a point x'^ where the ferromagnetic phase is preferred). The open symbols show 
that \E' predicts a transition at x'^ ~ 1.4, while the solid symbols, corresponding to \E's, 
approach each other but do not cross yet in the range of densities considered in the flgure. 
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FIG. 3. (Color online) Energy per particle of the polarized (grey triangles) and unpolarized (red 
diamonds) system obtained from the HNC/EL-optimized wave function of Eq. ([9]) including backs- 
flow correlations, in the framework of the FHNC theory. The solid and dotted lines are a guide to 
the eye. 



In any case, it is relevant to realize that the use of a Hamiltonian allowing the interaction in 
all three spin channels (up-up, down-down and up-down), together with a rich variational 
wave function, bring the transition point back from x'^ ~ 1 to higher values closer to the 
prediction of the Stoner model, x'^ ~ 1.66. 

The precise density at which our fully correlated model predicts the transition to the 
ferromagnetic state can be determined by extending the range of x values considered. In 
Fig. [3] we show our predictions for E^ /N and E''^'^ /N obtained from the optimized wave 
function \1/b of Eq. ([9]) including backflow correlations. As it can be seen, there is a wide 
range of densities where E^ /N and E^'^ /N are fairly close to each other, and thus estab- 
lishing a clear prediction of the transition point is delicate. Still, we clearly see that for 
X > 1.8 the polarized gas has lower energy than the unpolarized one, showing that at this 
point (and higher densities) the ferromagnetic phase is preferred. To the present accuracy 
we find that both curves cross at x'^ ~ 1.8, which establishes our prediction of the point 
where the ferromagnetic transition takes place for the hard-sphere Fermi gas. This transition 
occurs at a density that is close to the estimated freezing density of quantum hard spheres 
(xf ^ 1.95) [34] but slightly below. Our result for x'^ turns out to be siirprisingly close to 



a. 



the original predictions given in Ref. [2[, although recent experiments |5| with Li atoms 

11 



seem to indicate that in nature a pairing instability breaks the ultracold nature of the gas 
before the transition takes place. In any case and as stated above, one should recall that 
the mechanisms explored here are related but not equal to those found in the experiments 
which exploit the physical properties of the system close to a resonance. 

IV. SUMMARY AND CONCLUSSIONS 

To summarize, in this work we have analyzed the possible onset of a phase transition to 
a ferromagnetic phase for a two-component Fermi gas interacting through a repulsive hard 
sphere potential. Motivated by recent experimental and theoretical works, we have compared 
our results with existing predictions |6|, |l3| for a simplified Hamiltonian where only correla- 
tions between different spins are allowed. We have shown that, when properly renormalized 
to correctly account for p-wave scattering processes, already existing low-density expansions 
accurately reproduce Monte Carlo simulations of this simplified model. Based on that, we 
have shown that p-wave terms contribute significantly to the total energy per particle of the 
polarized and unpolarized systems at kpa ~ 1, and that therefore the simplified Hamiltonian 
for a fully repulsive interaction can not be used in this range. Finally, we have also presented 
calculations in the framework of the FHNC equations for a wave function including backflow 
correlations and an optimized Jastrow factor, showing that a ferromagnetic transition takes 



a 



place at kpa ~ 1.8, surprisingly close to the original prediction of Ref. J2|. 
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